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Abstract 

This contribution belongs to a combinatorial approach to hyperbolic geometry and 
it is aimed at possible applications to computer simulations. 

It is based on the splitting method which was introduced by the author and which is 
reminded in the second section of the paper. Then we sketchily remind the application 
to the classical case of the pentagrid, i.e. the tiling of the hyperbolic plane which is 
generated by reflections of the regular rectangular pentagon in its sides and, recursively, 
of its images in their sides. From this application, we derived a system of coordinates 
to locate the tiles, allowing an implementation of cellular automata. 

At the software level, cells exchange messages thanks to a new representation which 
improves the speed of contacts between cells. In the new setting, communications are 
exchanged along actual geodesies and the contribution of the cellular automaton is also 
linear in the coordinates of the cells. 

INTRODUCTION 

Hyperbolic geometry becomes more and more attractive, probably due to the 
irresistible aesthetic impression given by tihngs which are obtained in the hyper- 
bolic plane. An other reason to this attraction seems to be the fact that difficult 
problems can be solved more easily in the hyperbolic plane. This was first dis- 
covered in |22| . But an actual implementation in cellular automata was given in 
| 18I19| . This is the reason why the hyperbolic plane could be seen as a possi- 
ble tool for application purposes, see |14I2| . which are based on the approached 
defined in |^ and later generalized in 

First, in section 1, we remind a few features about the hyperbolic plane, 
in order the paper be self-contained. Then, in the section 2, we present the 
splitting method, which we introduced in |11I12| . with its tightly connected 
notion of combinatoric tilings. It introduces a general way to define a system 
of coordinates of the tiles. In section 3, we briefly remind the classical cases of 
the pentagrid and of the heptagrid. At this occasion, we remind the system of 
coordinates of the tiles for these tilings. Later, these coordinate systems will 
be considered as a hardware property of the implementation. In section 4, we 
remind the results of which allow to go continuously from the pentagrid to 
the heptagrid and conversely. In section 5, we introduce the new representation 
at the software level. The new system is closer based on the tree structure of 



the tiling. We describe a protocol for exchanging messages between cells from 
which we derive the properties which we announce: the connection between cells 
through geodesic paths and their computation is in linear or square time with 
respect to the length of the path. 

1. THE HYPERBOLIC PLANE 

We refer the reader to j24j for introductory material on hyperbolic geometry. In 
order to fix things and to spare space for this paper, we shall use Poincare's disc 
model of the hyperbolic plane, see figure 1, below. 




Figure 1 The Poincare's disc as a model of the hyperbolic plane. Here, p and q are 
parallel to £ and m does not cut £. 

The open unit disc U of the Euclidean plane constitutes the points of the hyper- 
bolic plane IH^. The border of U, dU, is called the set of points at infinity. 
Lines are the trace in U of its diameters or the trace in U of circles which are 
orthogonal to dU. The model has a very remarkable property, which it shares 
with the half-plane model: hyperbolic angles between lines are the Euclidean an- 
gles between the corresponding circles. The model is easily generalised to higher 
dimension, see |21I16| . 

This model gives a concrete realisation of hyperbolic geometry within math- 
ematics. It gives an example of a geometry where all the axioms of the Euclidean 
geometry are satisfied, except the parallel axiom. Here, this axiom is no more 
true: for any point A out of a line £, there are exactly two lines which pass through 
A and which are parallel to £. There are also infinitely many of them which pass 
through A and which do not intersect i, we call such lines non-secant. See 
figure ^ above, where both situations are illustrated. 

A direct consequence of this new property about parallel lines is that the 
sum of the angles of any triangle is less than tt. Indeed, the difference of this 
sum to TT is a measure: it is, by definition, the area of the triangle. 



Wc finish with this point by a last consequence: there is no notion of similarity 
in hyperbolic geometry. In particular, two triangles with the same angles are 
equal in this geometry, by contrast with what happens in Euclidean geometry. 

Recall that in the theory of tilings, a particular case consists in considering 
a polygon S and then, all the reflections of 5* in its sides and, recursively, of its 
images in their sides. All the polygons which are obtained in this way are said 
to be generated by S. We say that the set of polygons being generated by S is 
a tiling if and only if the interiors of two distinct polygons either coincide or 
do not intersect and, if any point of the plane is in the closure of at least one 
polygon of the set. In this case, we say that the tiling is generated by S. A tiling 
which is obtained in this way is also called generated by tessellation from S. 
We shall also say that the tiling is the grid defined by S. 

2. THE SPLITTING METHOD 

2.1. Basis of splitting and spanning tree 

Definition 1 — Consider finitely many sets Sq, . . ., Sk of some geometric 
metric space X which are supposed to be closed with non-empty interior, un- 
bounded and simply connected. Consider also finitely many closed simply con- 
nected bounded sets Pi, . . ., Ph with h < k. Say that the Si 's and Pi 's constitute 
a basis of splitting if and only if: 

(i) X splits into finitely many copies of Sq, 

(a) any Si splits into one copy of some P( and finitely many copies of Sj 's, 

where copy means an isometric image, and where, in condition (ii), the copies 
may be of different Sj 's. Si being possibly included. 

As usual, it is assumed that the interiors of the copies of Pi 's and the copies 
of the Sj 's are pairwise disjoint. 

The set Sq is called the head of the basis and the Pi 's are called the gener- 
ating tiles. 

Consider a basis of splitting of X, if any. We recursively define a tree A which 
is associated with the basis as follows. First, we split 5*0 according to the condi- 
tion {ii) of definition n This gives us a copy of say Pq which we call the root 
of A and which we call also the leading tile of Sg. In the same way, by the 
condition (ii) of definition ^ the splitting of each Si provides us with a copy 
of some Pi which we call the leading tile of Si. The splitting provides us also 
with ki regions, , . . ., Si^, which enter the splitting of Si. The regions which 
enter the splitting of 5*0 according to condition (ii) of definition ^ are called the 
regions of the generation 1. Assume that we have all the regions of generation n: 
Sn^, . . ., S„^^^. By definition, their leading tiles constitute the nodes of the gen- 
eration n. We split again these ^j's according to the condition {ii). We obtain 
rUn tiles which are called the tiles of the generation n-\-l and, for each S'„^ which 
is some Si, we have a splitting which is the isometric image of the splitting of Si, 
as above indicated. We say that the leading tiles of these copies of the splitting 



of Si are called the sons of the leading tile of . By definition, the sons of the 
leading tile of Sn^ belong to the generation n+1. The union of all the sons of 
the nodes of the generation n constitutes the nodes of the generation 

This recursive process generates an infinite tree with finite branching. This 
tree, A, is called the spanning tree of the splitting, where the splitting refers 
to the basis of sphtting So, ■ ■ ■ ,Sk- 

2.2. Combinatoric Tilings 

We come now to the following general definition: 

Definition 2 — Say that a tiling of X is combinatoric if it has a basis of 
splitting and if the spanning tree of the splitting yields exactly the restriction of 
the tiling to So, where So is the head of the basis. 

From |llll2j . we know that when a tiling is combinatoric, there is a poly- 
nomial which is attached to the spanning tree of the splitting. We have the 
following result: 

Theorem 1 — (Margenstern, |11I12| ') Let T be a combinatoric tiling, and de- 
note a basis of splitting for T by So, ■ ■ - , Sk with Po, . . . , Ph as its generating 
tiles. Let A be the spanning tree of the splitting. Let M be the square matrix 
with coefficients rriij such that rriij is the number of copies of Sj-i which enter 
the splitting of Si-i in the condition {ii) of the definition of a basis of splitting. 
Then the number of nodes of A belonging to the generation n is given by the sum 
of the coefficients of the first row of M" . More generally, the number of nodes 
of the generation n in the tree which is constructed as A but which is rooted in 
a node associated to Si is the sum of the coefficients of the row i+1 of AP^ . 

This matrix is called the matrix of the splitting and we call polynomial of 

the splitting the characteristic polynomial of this matrix, possibly divided by 
the greatest power of X which it contains as a factor. Denote the polynomial by 
P. From P, we easily infer the induction equation which allows us to compute 
very easily the number u„ of nodes of A on its level n. This gives us also the 
number of nodes of each kind at this level by the coefficients of M" on the first 
row; we use the same equation with different initial values. Sequence {u„}„giv 
is called the recurrent sequence of the splitting. 

Now, as in |9|1 1 12 , number the nodes of A level by level, starting from the 
root and, on each level, from the left to the right. Consider the recurrent sequence 
of the splitting, {m„}„>i: it is generated by the polynomial of the splitting. As 
we shall see, it turns out that the polynomial has a greatest real root (3 with 
(3 > \. Sequence {un}n>i is increasing and as well known, see |4I7I for instance, 

k 

it is possible to represent any positive number n in the form n = Oi.Ui, where 

Oj € {0..[/3J }. The string a/c . . . ao is called a representation of n. In general, the 
representation is not unique and it is made unique by an additional condition: 
we take the representation which is maximal with respect to the lexicographic 



order on the words on {0..6} where b — [/3J . The set of these representations is 
called the language of the splitting. 

In this paper, we consider only the case when we have a single generating 
tile, i.e. when h = 1. 

A lot of tilings are combinatoric, see references in |lUll3j . not only tilings in 
the hyperbolic plane but also in Euclidean spaces. Moreover, in several cases of 
tilings of the hyperbolic plane, it turns out that the language of the splitting is 
regular. 

In the next section, we exemplify the application of the method on the pen- 
tagrid, the tiling of IH^ which is generated by the rectangular regular pentagon. 

3. THE CLASSICAL CASE OF THE PENTAGRID 

Here, the basis of the splitting consists of two regions, Q and R^. The first region, 
Q, is a quarter of the hyperbolic plane, the south-western quarter in figure 2. 
It constitutes the head of the splitting. The second region, i?3, which we call a 
strip, appears in figure 2 as the complement in Q of the union of the leading 
tile Pq of Q together with two copies of Q, Ri and i?2, see the figure. Region i?i 
is the image of the shift along side 1 of Pq which transforms a vertex of Pq into 
the other vertex of the side. Region R2 is obtained from Q by the shift along 
side 4 of Pq- 
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Figure 2 The splitting of Q which is associated to the pentagrid. Notice the construc- 
tion of the spanning tree. 

This gives us the splitting of Q. We have now to define the splitting of R3. 

Its leading tile is provided us by the reflection of Pq in its side 4, say Pi. Call 
R[ the image of Q by the shift along the side 5. It is also the reflection of Pi in 
the diagonal of Q. The shift along the side li of Pi transforms P'l into 5*1. Now, 
it is not difficult to see that the complement 5*2 in P3 of the union of Pi and 5*1 



is the image of i?3 under the shift along the side 5 of Pq- And so, R3 also can 
be split according to the definition. 

The exact proof that the restriction of the tiling is in bijection with the 
spanning tree of the splitting is given in |18I20| . 

In the case of the pentagrid, the number of nodes of the spanning tree which 
belong to generation n is /2n+i where {fn}n>i is the Fibonacci sequence with 
fi — I and /2 = 2. This is why the spanning tree of the pentagrid is called the 
Fibonacci tree, starting from |18I2(J| . 

From figure 2, and arguing by induction, it is plain that the Fibonacci tree 
is constructed from the root by the following two rules: 

- a 3-node has three sons: to the left, a 2-node and, in the middle and to 
the right, in both cases, 3-nodes; 

- a 2-node has 2 sons: to the left a 2-node, to the right a 3-node. 

and starting from the axiom which tells us that the root is a 3-node. 

Starting from 9 , we call this tree the standard Fibonacci tree, because there 
are a lot of other Fibonacci trees. As it is proved in there are continuously 
many of them, and each one is associated to a different splitting. We refer the 
interested reader to [5] or to our previous surveys |1UI13) . 

In PI, a new and more efficient way is defined to locate the cells which lie 
in the quarter, which is the prefiguration of the algebraic side of the splitting 
method later introduced in 12 . We number the nodes of the tree with the help 
of the positive numbers: we attach 1 to the root and then, the following numbers 
to its sons, going on on each level from the left to the right and from one level 
to the next one, see figure 3, below. 
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Figure 3 The standard Fibonacci tree: 

above a node: its number; below: its standard representation. 

Notice that the first node of a level has a Fibonacci number with odd index as its 
number. 



Figure 4 Above, the ternary heptagrid by one central tile and seven copies of So- 
Below, the pentagrid by one central tile and five copies of So ■ 

This numbering is the basis of the system of coordinates which was devised 
for representing the tiles of the pentagrid starting from As an example, this 
system was used in implementations as in . The system is the following: as 



splits into four copies of 5o, it is enough to attach a number in {0..3} to 
each copy and then, inside a copy, to use the just described numbering. 

Now, we know that positive integers can be represented with the help of 
the Fibonacci sequence, again see \7\ for instance, which is a very particular 
case of what was indicated in section 2. Indeed, any positive integer is a sum of 
distinct Fibonacci numbers. The representation is not unique, in general. It is 
made unique by an additional condition: we take the longest representation or, 
equivalently, we do not accept representations which contain the pattern '11'. 
The longest representation is called the standard Fibonacci representation. The 
language of these representations is regular, which is a classical result in lan- 
guage theory, see [7], and we call coordinate of a node the standard Fibonacci 
representation of its number in A. 

4. THE FIBONACCI CARPETS 

As noted in |T7) . it is natural to wonder whether it is possible to transport a 
configuration within the pentagrid into the heptagrid and conversely, without 
big distortions. In |17| . we showed that this is possible and, in the quoted paper, 
we indicate a solution which is based on the notion of carpet. 

The idea can be summarized as follows. 

In the hyperbolic plane, let us fix an angle a with a < tt. Then the angular 
sectors of angle a always contain a half-plane. As it is easy to see that M^is a 
growing sequence of half-planes, i. e. each half-plane is contained in its successor 
in the sequence, we can expect that has the same property with the angular 
sectors. And this is the case. 

From this property, it is not difficult to conclude that we can do this with 

TT 

quarters: it is the case when a = — . As proved in J7], it is not difficult to con- 
struct a growing sequence of quarters {i?„}„e^ possessing the following prop- 
erty. Denoting the Fibonacci tree rooted in the leading tile of Hn by F„, we 
require that the root of Fn is the middle son of the root of F„+i. Of course, we 
can do the same in the ternary heptagrid. This time there is a growing sequence 
of angular sectors {Kn}ni^z: which covers IH^ . These sectors have the angle /3 
between the mid-point lines which define a sector in the ternary heptagrid, see 
| 2I17 | . Here also, as proved in ^7] , denoting the Fibonacci tree associated to Kn 
by G„, we have that the root of G„ is the middle son of the root of G„+i. 

Now, the bijection between Fn and G„ is easy: it is enough to take tiles with 
the same coordinate with respect to the corresponding roots. 

The carpet representation is defined as follows. 

Consider a tile T of the pentagrid. There is a unique n such that T G Fn 
and T ^ Fn-i- This n is the first coordinate of T. The second coordinate of T 
is its coordinate in F„. Denoting the coordinates of T by (n, v), this define an 
injection from the tiles of the pentagrid into ZxlN . We can define coordinates 
in the ternary heptagrid in exactly the same way. The expected bijection consists 
in putting in correspondence tiles with the same coordinates. As proved in |17| . 



this defines a continuous bijection from the pentagrid into the heptagrid whose 
inverse is also continuous. 

In order to give their full strength to the above results and to enlighten the 
next section, we define the notion of distance in our context. 

A path joining the tile Tq to the tile Ti is a sequence {Si}o<i<n of tiles such 
that: 

(i) So — Tq, Sn = Ti, 

(a) for each i, with < i < n, Ti and T^+i have a common side; the 
number n is called the length of the path. 
We call distance between Tq and Ti, denoted by dist(To,Ti), the shortest 
length for a path joining Tq to Ti. It is very easy to show that the distance which 
we have just defined satisfies the traditional axioms of a metric. A geodesic path 
from To to Ti is a path joining Tq to Ti whose length is dist(To, Ti). 
We have the following property: 

Lemma 1 In a Fibonacci tree, there is a single geodesic path from the root to 
another cell: it is the branch issued from the root which goes through this cell. 

Proof. It is an easy induction on the levels of the tree: the level n is the set of 
nodes of the tree which arc at distance n to the root. | 

5. COMMUNICATIONS BETWEEN THE CELLS 

In order to simultaneously deal with both the pentagrid and the heptagrid, we 
denote by p the number of sides of the polygon which generates the considered 
tessellation. Accordingly, p = 5 for the pentagrid and p — 7 for the heptagrid. 

In this section, the tiles are called cells, as we assume that they have a finite 
automaton, the same for all the tiles, at their disposal. When we shall consider 
a cell within a tree and from the point of view of the tree connections, we shall 
use the term node to denote the tile of the cell. As usual for cellular automata, 
each cell is in contact with its neighbours. Here, the neighbours of a cell T are 
the cells whose tile shares a side with the tile of T. 

5.1. The problem 

Consider the following problem which generalizes the situation considered in 
|17|. A cell c sends a message to all cells in order to establish communication 
with the cells which are interested by the message. Such a message will be called 
public. When a cell answers to the message of another which is known, then the 
message is private. As we shall see, the mechanism of routing the messages is 
not the same depending on whether the conveyed message is public or private. 
In this paper, we consider only the part of the transmission which is concerned 
by the computation of the address of the cell to which a message is sent. We do 
not take into account the size of the message: this raises the question how the 
message is produced, what are the storage capabilities of the cell and so on. As 
the resources of a cell of a cellular automaton are bounded, the transmission of 
messages must be dealt with by the hardware. We shall consider that the cell 



only indicates whether it has a message or not. In the case of a message, the cell 
also indicates whether it is public or private. In this context, the only complex 
thing is the computation of the addresses. 

In order to be contacted for a reply, the sender c must send its address to every 
cell. This will allow a receiver d of the message to send back an acknowledgment 
or a message accepting the contact to c only. As the address for distant cells 
may exceed the resource of a cell, the transmission of addresses is also performed 
by the hardware. However, the cells can contribute to the computation of the 
address. This is what we assume and we deal with this problem. 

In the systems of coordinates of |21 or ^Tj, the addresses being absolute 
entail that the connection will not be along the shortest path. As proved in 
and in '17", the path has a length which is linear in the addresses of c and d. In 
the solution which we now suggest in this paper, the path will still be linear in 
the addresses of c and d, but now, it will also be a geodesic path between the 
cells. 

5.2. A solution 

The idea is to replace absolute coordinates by relative ones and to obtain the 
address of d in the shortest time. 

There is a fixed system of coordinates, for instance, the system described in 
12]. The system is based on a central cell Tq with coordinate (0, 0) and p angular 
sectors around To, covering IH^. The p Fibonacci trees associated to the angular 
sectors together with Tq exactly contain all the tiles of the tessellation. 

For each cell T with T Tq, its coordinate is (/, i^) where / G is the 

number of the Fibonacci tree to which it belongs and ly is the coordinate of T 
in this tree. This system is used by the hardware in order to locate the cells. 

Now, when the cell c emits a public message, c considers itself as the central 
cell of a system of coordinates. Also, in order to obtain a fast computation, the 
message arrives to each cell e together with the address of e relative to c. Our 
problem will be to manage the updating of the computed address, each time 
the message arrives to a new cell. The goal is to realize this operation as fast as 
possible. 

We solve the problem by introducing a new system of coordinates which will 
be used both as an absolute system and as a relative one. We first describe the 
system, which we call the software system as it is used by the cells during 
their computations. Then we describe the protocol of the message management 
and the computations it entails for the cells. 

The software system 

We introduce new absolute coordinates for the location of a cell. We describe 
it in a way which is common for the penta- and the heptagrid. A bit further, we 
shall indicate another solution which is specific to the pentagrid. 

These coordinates consist of digits in [l..p], with p e {5, 7}, each one repre- 
senting an arc from a node to one of its neighbours. Except for the central cell 
which has no father, all the other cells have a father, and the arc joining a node 
to its father has the number 1. The other neighbours of the corresponding cell 



are joined by the arcs numbered from 2 up to p, while counter-clockwise running 
around the node. This also holds for the roots of the Fibonacci trees which are 
considered as the sons of the central cell. From this principle, we remark that 
a given arc does not receive the same number, depending on which sense it is 
crossed. For instance, an arc 1 goes from a node to its father. For the reverse 
direction, the number of the same arc is 2, 3 or 4, for the pentagrid, 3, 4 or 5, 
for the heptagrid, depending on the statuses of the nodes which are at the ends 
of the arc. In order to avoid problems connected with this difference in the de- 
notation of an arc, we merge both notations in a unique one as follows. Consider 
an arc from the side 7 of a cell c, in the case of the heptagrid. On the other side 
of this edge, the same side is numbered 2 or 3, depending on the status of the 
other cell d which shares the side with c. And so, if we go from c to d, the side 
will be denoted by (7, 2) or (7, 3). If we go from d to c, the arc will be denoted 
by (2, 7) or (3, 7) respectively. We introduce a mirror operation on such couples: 
let S = [a, (3) with a,l3 G {l..p}. We denote the mirror oi S by S = {P,a). In 
d = (a,/3), we say that a is the input and (3 is the output. We define corre- 
sponding operators on 6 by: a = {S)i and /3 — {d)o- Now, the definition of the 
mirror will allow us to shorten the notations: if we consider all possible couples of 
input and output in this order, the set of couples is divided into two equal classes 
which are the mirror of each other. Accordingly, we introduce the notations of 
tabled 

For the digit which corresponds to the connection from the central cell to a 
neighbour or for the opposite connection, we have a similar convention. However, 
the connections are of the form (i, 1) with i G {l.-p}. To avoid confusions, we 
encode these couple with bold digits in {l..p}. 
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The central cell has as its coordinate. 

For the other nodes, the first digit indicates in which Fibonacci tree the 
node is. We call this first digit the sector. Next, the other digits are established 
according to the indicated principles. The considered sequence can be seen as 
a variant of the system of digits which occur in the linear algorithm defined in 

The protocol 

The same system is used for the exchange of messages. 



When a cell c sends a public message, it considers itself as the centre of 
a similar system of coordinates which we call relative by contrast with the 
absolute system. The relative address of c is 0. The message is sent to all the 
other cells of the plane. The other cells behave as relays: they receive the message 
sent by c together with a sequence of digits So -Si which is the sequence of arcs 
going from c to the relay; the relay forwards the message to its sons, within the 
relative tree rooted in c, and transmits them the sequence So..6iSi^i. Of course, 
for each son of the relay, J^+i is the arc joining the relay to the son. Accordingly, 
the work of the relay is only to append But here, we require more: the i5i's 
must denote arcs in the absolute system of coordinates. Now, the computation of 
Si-^-l, as an absolute arc, is easy. Assume that Si — (ao,/3o) and Si+i = (ai,/3i). 
We may assume that each cell knows what is the number of its edges in the other 
cell of the tiling which shares the same side. This information can be contained 
in a table output which is memorized in the states of the cell. Then, denoting a 
son by s, with s € {0..2}, from the leftmost to the rightmost son, ai and /3i are 
given by: 

ai = l + ((/3o-l) + 2 + ^^ + 5- st^)modp, 

(1) 

/3i = output{ai). 

where str is the status of the relay in the relative tree. 

Now, consider the case of a private message. It can always be considered as 
a reply to a message, either public or private. 

If a relay d wishes to reply to the message, it knows its own address from 
the point of view of c. But the elements of the address are absolute arcs. And 
so, backward taking the same path leads to c: the structure of the digits allows 
to perform this very easily. 

Indeed, for the reply, d will ask the hardware to copy the address sequence 
received from the previous relay which is the address sent by c. The cell d knows 
that the first relay to which it has to send the reply is on the other side of its 
edge And so, as d used the information contained in Si it asks the hardware 
to prepare two sequences to be transmitted to c: So--Si-i and Si. The relays will 
interpret these sequences as stacks whose top is on the right-hand end of the 
sequence, considered as a word. Accordingly, when a relay receive the reply with 
the following two sequences: So..Sj and it asks the hardware to pop the 

top of the first sequence and to push it on the second sequence. And so, it will 
send the following two sequences So--Sj-i and 5j..5i to the cell which is on the 
other side of its edge {Sj)o- 

It is easy to see that possible further exchanges between c and d will use the 
same mechanism. Both the strings will be sent with exchanging the roles each 
time the receiver is reached. In the direction from c to d, the relays read {S)i, 
where Si is the top of the first sequence, while in the direction from d to c, the 
relays read {S)o, where Si is also the top of the first sequence. Of course, there is 
a bit also sent with the sequences which indicates which part of the information 
should be read. 



An example 

First, consider the example in the heptagrid, see figure|Sl 




Figure 5 Illustration of the example in the heptagrid. The digits in light hue are in 
relative coordinates. 




Figure 6 Illustration of the example in the pentagrid. 



The sender is the cell with coordinate 312332. Consider the cell 2331332. It is 
not difficult to see that the address of d with respect to c is 23362332, see figureEl 



while it is 13232332 in relative coordinates with c as the central cell. The reverse 
path is 23326332 in absolute arcs: it is the reverse of the path from c to d. 

Note, that if we take the same example in the pentagrid, the coordinates are 
a bit changed. The absolute coordinate of c is 312332, the absolute coordinate 
of d is 2331332, but now, the relative coordinate of d from c is 13232332, and 
the same path in absolute coordinates is 33232332. 

Note that when a cell is inside a tree defined by the central cell, the digits of 
its coordinates are in {1..3}. The occurrence of other coordinates and especially 
of over-lined digits indicate either the crossing to another tree or sub-tree or a 
backward segment of the path. 

On the example of figure El note that the path from c to d is not exactly 
the same if we compare this path when in the pentagrid, when in the heptagrid. 
The crossing from one tree to the neighbouring one does not occur in the same 
way. This can be explained by the fact that in the heptagrid, each cell has more 
neighbours than in the case of the pentagrid and, sometimes, this gives a shorter 
path as we can see by comparing the figures[Sland|Sl If we consider the central cell 
of these figures as an ordinary node, we can see that the phenomenon indicated 
by the figures also occurs inside trees. 

A special system for the pentagrid 

TT 

In the pentagrid, due to the interior angle of — , we can use the following 

property, first noted in [201 and also more especially studied in [131. We can 
define a tihng d la Wang in the pentagrid by giving numbers in {1..5} to the 
edges of the pentagons by applying the following process. 

We start from a given pentagon on which we fix the numbering: we choose 
at random one of its edges which receives the number 1 . Then, we number the 
other sides by increasing numbers when running around the pentagon, counter- 
clockwise. Then this numbering fixes the numbers of all other pentagons in such 
a way that all edges are always given the same number in the pentagons to which 

they belong. As this property is not true for the ternary heptagrid due to the 
27r 

angle — , the construction works for the pentagrid only, 
o 

We define the absolute and relative coordinates as previously with this dif- 
ference that the digits are now in {1..5}, as edges have the same number both 
from the input and the output sides. However, the algorithm to compute the 
digits which are appended is a bit different. We have to take into account that 
the numbering has the following property: if in a pentagon c the number are in- 
creasing when run counter-clockwise, they are decreasing for all the neighbours 
of c when the running is still performed counter-clockwise. Accordingly, there 
is a notion of orientation. By definition it is positive for the central cell and 
it is changed in the opposite orientation each time an edge is crossed. In this 
new setting, we assume that the cells know their orientation. In this case, if So 
is the top of the sequence, if sU is the status of the present node in the relative 
tree and if or S {— 1,-M} is the absolute orientation of the cell, the new digit Si 



corresponding to the son s in the relative tree, with s G {0..2}, is given by the 
following formula: 

6i = l + ((So-1) + orx{2+s-str))mod 5. (2) 

If we consider the example of figure El the cell c has the coordinate 324142, 
the cell d has the coordinate 2421413. The path from c to c? in absolute coordinate 
is 242131413. 

Theorem 2 The algorithm of section 5 using table ^ and formulas ^ allows 
to establish geodesic communications between any cells of the pentagrid or the 
ternary heptagrid. The similar algorithm based on the natural colouring of the 
pentagrid and formula^ does the same for the pentagrid. 

Proof. The 'geodesic' part of the proof is a direct consequence of lemma It 
comes from the fact that the address which arrives to any cell contacted by c is 
built along a branch of a Fibonacci tree. B 

5.3 Complexity issues 

Now, we have to look at the complexity part of the problem. First, we consider 
the software aspect: the complexity of the work of the cellular automaton. Then 
we look at the hardware side of the issue. 

Software side 

We have to look exactly on the conditions of the working of the messages. We 
remember that the role of the cellular automaton consists only on computing a 
new digit while forwarding a public message. 

Accordingly, each cell of the cellular automaton indicates to which neighbour 
a received message has to be forwarded. Each cell also has to compute the part 
of the transmitted address which has to be appended to the computed address 
in the case of a public message. In the case of a private message, the mechanism 
of popping a stack and pushing the popped information on the other stack can 
be handled by the hardware. This mechanism detaches the popped digit which 
is transmitted to the cell together with the bit which indicates the direction 
of reading. The correct input is obtained which allows the cell to get the right 
output. The table involved in the computation of formula^can be stored in the 
states of the automaton. 

In this case, the computation performed by each cell is constant at each top 
of the clock and so, the above algorithm is linear in time with respect to the 
distance between the two cells: 

Theorem 3 The execution of the algorithm considered in theorem^is linear in 
time in the distance between the two communicating cells. 

We also note that at each top of the clock, each cell receives at most p mes- 
sages as it has p neighbours. As p is a constant of the automaton, it may be 



contained in the states of the cells. Accordingly, we may assume that the com- 
putation of the new digit for the address of a public message is simultaneously 
performed for all the p messages received by the cell between two consecutive 
tops of the clock. 

However, we have to take into account that a cell may receive several messages 
for the same destination. It is not difficult to consider cases when this number 
exceeds the resource of a cell: imagine that all the cells at distance k reply to c. 
If all messages to the same cell are freely transmitted, there are pf2k+i messages 
arriving at c, where /„ is the n*^^ term in the Fibonacci sequences defined by 
/o = /i = 1. For this reason, we require that when more than one message arrive 
at a cell d for the same output, then only one message is transmitted. The others 
are stored at d in an appropriate stack provided by the hardware. At each top 
of the clock, the cell gives priority to the messages waiting in the stack which 
works in a FIFO mode. 

In this way, messages are guaranteed to reach their goal and there cannot be 
dead-locks in the system. This includes the possibility for any cell cell relaying 
private or public messages may also decide to send a public message. But, there 
may be delays, possibly important, during the transmission. However, if we as- 
sume that communications with distant cells are rare, the above situation will 
not occur and delays will also be rare. 

Indeed, assume that the probability for two cells c and d to start a com- 
munication is (7e~^'*"'*(°''^), where C and A are appropriate constants. Let Pk 
be the probability for c to have a communication with a cell at distance k. 
Then, as the number of cells at distance k from c is pf2k+i, we have that 
Pk = Cp/2fe+ie--^<^'"*(=''^). Now, we know that /afe+i ~ Cl(/?)^ where /? is the 
dominant root of X"^ — 3X -|- 1 and C\ is an appropriate constant. A simple com- 

3 ^- \/5 

putation gives us that fi = — ^ — < e. Accordingly, Pk < Dj'', with D = CC\p 

and 7 = If we denote by P<k and P>ft the probabilities for c to have a com- 
munication with a cell at a distance at most k or greater than k respectively, we 

have that P>fc < D- and, consequently, P<k > 1 — -D- . We fix D and 

1 — 7 - 1 — 7 

oo 

A so that ^Pfe = 1. This computation shows us that our hypothesis of rare 

fc=0 

distant communications is sound. Accordingly, we may consider that the size of 

the delays and their occurrence as bounded. 

At last, another important aspect of the considered system is that the same 
message sent by one cell arrives at a given cell exactly once. This comes from 
the bijection between the tiling and the union of p Fibonacci trees together with 
the central cell. 

Hardware side 



At the hardware level, we note that going from one cell to another involves 
the computation of the new absolute address thanks to the given absolute ad- 



dress. With the system of coordinates which is indicated at the beginning of 
section 5.2, taken from |9I14| . the just considered computation is hnear with 
respect to the given coordinate. Now, from going to one end of a path to the 
other, the computation is then quadratic as the length of the path is hncar in 
coordinates of the ends of the path. 

Now, consider the mechanisms which we introduced in section 5.2. We noted 
that the management of the stacks involved in these mechanisms must be per- 
formed by the hardware. As the total length of the stack is always equal to the 
length of the path, for each relay, the management of the stacks is linear in the 
length of the path. And so, the cost of this management is quadratic in the 
length of the path. And so, at the hardware level, the new system has exactly 
the same complexity as the absolute system based on the coordinates defined in 
P]. As we may assume a bigger speed for the hardware, the quadratic part of 
the hardware manipulation may be not higher than the linear cost at the level of 
the cellular automaton for distances which are not higher than the ratio between 
software and hardware speeds. 

For what is messages, the transmission from a cell c to a neighbour consists 
in copying the message present in c to the appropriate support of d. If Af is the 
length of the message, the cost of the copy is Mp, where p is the time cost of an 
elementary hardware operation. Accordingly, if AI is the length of a message, the 
complexity of the transmission from a cell c to a cell d is Mpxdist(c, d). And so, 
the estimate of the total cost of conveying a message depends on the comparison 
between M and the distance run by the message. If the distance is big, then the 
computation of the addresses is more expensive than the pure transmission of 
the message which is copied from one cell to the next relay. If the length of the 
message is bigger than the address, then this length is the main factor in the 
complexity estimation. In fact, the complexity is given by the following result: 

Theorem 4 The overall complexity of transmitting a message of length M with- 
out delay from a cell c to a cell d is dist(c, d) (l + pA/+pdist(c, d)) , where p is the 
time cost of an elementary hardware operation, the time cost of an elementary 
software operation being 1. 



CONCLUSION 

As a conclusion of this paper, we hope that the results of this paper can be 
of help for implement purposes. 

It is time to refer here to the few applications of the results quoted from 
jyill2| . These works gave way to one application in human-machine interface, 
see P], where a colour palette is presented to allow a user to find easily hues 
starting from fundamental colours. 

We think that there are possible other applications. As indicated in 14 , one 
of them could be the representation of the Internet. Another possible application, 
also mentioned in H^l is special relativity theory. Section 5 could be of interest 
in both directions. 



References 



1. Chelghoum, K., Margenstern, M., Martin, B., Pecci, I.: Locating Points in the Pen- 
tagonal Regular Tiling of the Hyperbolic Plane, proceedings of SCI'2003, Orlando, 
USA, July, 27-30, 2003. 

2. Chelghoum, K., Margenstern, M., Martin, B., Pecci, I.: Cellulax automata in the 
hyperbolic plane: proposal for a new environment, proceedings of ACRI'2004, 3305, 
Amsterdam, Ncthorland, October, 25-27, 2004. 

3. Chelghoum, K., Margenstern, M., Martin, B., Pecci, I.: Palette hyperbolique : un 
outil pour interagir avec des ensembles de donnees, proceedings of IHM'2004, 
Namur, Belgium, August, 30 - September, 3, 2004. 

4. Praenkel, A. S.: Systems of numerations, American Mathematical Monthly, 92, 
1985, 105-114. 

5. Harborth, H.: Concentric cycles in mosaic graphs, 3, 1990, 123-128. 

6. Herrmann, F., Margenstern, M., Morita, K.: Cellular automata in the hyperbolic 

plane: implementing a polynomial algorithm to solve SAT problem. Proceedings 
of CA 2000, 6"^ IFIP WG1.5 meeting, Osaka, Japan, August, 21-22, 2000, 

7. Lothaire, M.: Algebraic combinatorics on words, Cambridge University Press, Cam- 
bridge, 2002. 

8. Margenstern, M.: Cellular automata in the hyperbolic plane. Publications du 
GIFM, N° 99-103, Metz, France, 1999, ISBN 2-9511539-4-5. 

9. Margenstern, M.: New tools for cellular automata in the hyperbolic plane. Journal 
of Universal Computer Science, 6(12), 2000, 1226-1252. 

10. Margenstern, M.: Cellular automata in the hyperbolic plane. A survey, Romanian 
Journal of Information Science and Technology, 5(1-2,), 2002, 155-179. 

11. Margenstern, M.: A contribution of computer science to the combinatorial ap- 
proach to hyperbolic geometry, SCI'2002, 5, Orlando, USA, July, 14-19, 2002, 
ISBN 980-07-8150-1. 

12. Margenstern, M.: Revisiting Poineare's theorem with the splitting method, 
Bolyai'200, International Conference on Ceometry and Topolgy, Cluj University 
Press, Cluj-Napoca, Romania, October, 1-3, 2002, ISBN 973-610-205-X. 

13. Margenstern, M.: Cellular Automata and Combinatorie Tilings in Hyperbolic 
Spaces, a survey, proceedings of DMTCS'2Q03, 2731, Dijon, France, July, 7-12, 
2003. 

14. Margenstern, M.: Implementing Cellular Automata on the Triangular Grids of the 
Hyperbolic Plane for New Simulation Tools, proceedings of ASTC'2003, Orlando, 
USA, March, 29 - April, 4, 2003. 

15. Margenstern, M.: Tiling the Hyperbolic Plane with a Single Pentagonal Tile, Jour- 
nal of Universal Computer Science, 8(2). 2002, 297 316. 

16. Margenstern, M.: The tiling of the hyperbolic 4D space by the 120-cell is combi- 
natorie, Journal of Universal Computer Science, 10(9), 2004, 1212-1238. 

17. Margenstern, M.,: A New Way to Implement Cellular Automata on the Penta- and 
Heptagrids. Journal of Cellular Automata, 1, 2006, 1-24. 

18. Margenstern, M., Morita, K.: NP problems are tractable in the space of cellular 
automata in the hyperbolic plane. Publications du LIFM, N° 98-101, Metz, Prance, 
1998, ISBN 2-9511539-3-7. 

19. Margenstern M., Morita K., A Polynomial Solution for 3-SAT in the Space of 
Cellular Automata in the Hyperbolic Plane, Journal of Universal Computations 
and Systems, 5, (1999), 563-573. 



20. Margenstern, M., Morita, K.: NP problems are tractable in the space of cellular 
automata in the hyperbolic plane, Theoretical Computer Science, 259, 2001, 99- 
128. 

21. Margenstern, M., Skordev, G.: Tools for devising cellular automata in the hyper- 
bolic 3D space, Fundamenta Informaticae, 58(2), 2003, 369-398. 

22. D. Morgcnstcin and V. Kroinovich, Which algorithms are feasible and which are 
not depends on the geometry of space-time, Geocombinatorics, 4 3 (1995) 80 97. 

23. Martin, B.: VirHKey : a VIRtual Hyperbolic KEYboard with gesture interaction 
and visual feedback for mobile devices, accepted to MobileHCI'2005, 7th In- 
ternational Conference on Human Computer Interaction with Mobile Devices and 
Services, Salzburg, Austria, September, 19-22, 2005. 

24. Meschkowski, H.: Noneuclidean Geometry, Acad. Press., New- York, 1964, (trans- 
lated by A. Shenitzer). 



